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We present a theory of large-scale dynamo action in a turbulent flow that has stochastic, 
zero-mean fluctuations of the a parameter. Particularly interesting is the possibility 
of the growth of the mean magnetic field due to Moffatt drift, which is expected to 
be finite in a statistically anisotropic turbulence. We extend the Kraichnan-Moffatt 
model to explore effect s of finite memory of a fluctuations, in a spirit similar to that of 
Sridhar fc Singhl ( 20141) . hereafter SS14. Using the first-order smoothing approximation, 
we derive a linear integro-differential equation governing the dynamics of the large-scale 
magnetic field, which is non-perturbative in the a-correlation time T a . We recover earlier 
results in the exactly solvable white-noise (WN) limit where the Moffatt drift does not 
contribute to the dynamo growth/decay. To study finite memory effects, we reduce the 
integro-differential equation to a partial differential equation by assuming that the r a 
be small but nonzero and the large-scale magnetic field is slowly varying. We derive 
the dispersion relation and provide explicit expression for the growth rate as a function 
of four independent parameters. When r a ^ 0, we find that: (i) in the absence of the 
Moffatt drift, but with finite Kraichnan diffusivity, only strong a-fluctuations can enable 
a mean-field dynamo (this is qualitatively similar to the WN case); (ii) in the general 
case when also the Moffatt drift is nonzero, both, weak or strong a fluctuations, can 
lead to a large-scale dynamo; and (iii) there always exists a wavenumber ( k ) cutoff at 
some large k beyond which the growth rate turns negative, irrespective of weak or strong 
a fluctuations. Thus we show that a finite Moffatt drift can always facilitate large- 
scale dynamo action if sufficiently strong, even in case of weak a fluctuations, and the 
maximum growth occurs at intermediate wavenumbers. 
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1. Introduction 


The magnetic fields observed in various astrophysical bodies, such as the planets, 
the Sun, stars , galaxies etc, are believed to be self-sustained by turbulent dynamos 


dMoffatt 19781: Parker 1979: Krause & Radlei 198(]| Ruzmaif 

.in. Shukurov & Sokoloff 

Il988t Kulsrudl 12004: Brandenburg & Subramanianl 20051 Jonesl 

20 111). In an electrically 


conducting plasma, conversion of the kinetic energy into the magnetic energy, without any 
electric current at infinity, is known as the dynamo action, which leads to amplification 
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of weak seed magnetic field. Magnetic fields exhibit coherence over a range of scales, 
from smaller to much larger than the outer scale of the turbulence. Systems lacking 
ordered motion, such as the clusters of galaxies, predominantly host a fluctuation dynamo 
dMurgia et all 2004 Vogt & Enfilinl 200.4 Kucher fc Enfilinl 2011 : Bhat fc Subramaniaul 
121) 1 .'ill , whereas those w ith large-scale motion, such as t h e Sun, galaxies etc, also su pport 
a large-scale dynamo ( Charbonneau 20ld iBec aiml; IChamandv et ~all l2013allbl. I20l4 
and references therein). Dynamo origin of the galactic magnetic field seems unchallenged. 
Helical turbulence has been considered to be the key driver for large-scale dynamo, which 
could operate even in systems without mean motion. For example, both, analytical calcu¬ 
lations and numerical simulations, reveal that even in the absence of differential rotation 
or mean motion, large-scale m agnetic field grows due to helically f orced turbulence by 
what is known as an a 2 -effect ( Brandenburg fc Subramanian 120051 ). However, it is not 
clear whether astrophysical turbulence has a mean helicity that is large enough to sustain 
such a large-scale turbulent dynamo. This brings us to a following natural question: could 
the large-s cale magnetic fields grow if the helicity of the turbulence vanishes on average? 
Kra ichnan (: 197 6h was the first to study this problem by considering a (which is a measure 


of mean kinetic helicity of the turbulence) as a stochastic variable, with zero mean, and 
demonstrated that the a fluctuations lead to a decrement of the turbulent diffusivity, 
and if sufficiently strong, they coul d give r i se to the growth of mean magnetic field 
by the process of negative diffusion. iMoffattl (1978) generalized this model to include a 


statistical correlation between the fluctuating a and its spatial gradient, and found that 
this contributes a const ant drift vel o city t o th e dynamo actio n, but it does not affect the 
dynamo condition. Both Kraichnan! ( 19761) and lMoffattl ( 1978 ) have essentially considered 


white-noise a fluctuations, as has been explicitly shown by SS14 who extended previous 
studies by also studying the memory effects, giving rise to new interesting mechanisms. 

The SS14 model is limited to fairly low wavenumbers because their FOSA calculation of 
the mean EMF ignored turbulent resistivity. The present paper remedies this, obtaining 
a new expression for the mean EMF which predicts a significant modification of dynamo 
action at intermediate and large wavenumbers. We show that the inclusion of the resistive 
term in the present work to determine mean EMF is a nontrivial extension of SS14 model, 
leading to even qualitatively new predictions for the growth r ates at intermediate and 
large wavenumbers. Considering deterministic Roberts flows, iRheinhardt et al ( 2014 ) 


discussed the role of finite memory for the existence of the dy namo. Based on numerical 
experi ments of passive scalar diffusion and kinematic dynamos. iHubbard fc Brandenburg 


(2009) investigate turbulent transport where the turbulence possesses memory. These 


works highlight the importance of the memory effects on the turbulent transport processes 
and demonstrate that the turbulent transport coefficients can be significantly different 
from those cases where the correlation times of the turbulence are nearly vanishing. 

A number of previous studies have exploited the idea of a fluctuations to study 
the large-scale dynamo mechanism in wide variety of contexts. As many astrophysical 
sources possess differential rotation, the focus has been on th e understanding o f large- 


1997: Sokolov 19971: 

Silant’ev 

Richardson & Proctor 

2012ft. 


20001 : Proctor 2007 


In context to the solar dynamo 


20121 H ur fc Subramanianl 1200 
Silant’evl 


200 ' 


proposed a new dynamo mechanism by considering an inhomogeneous distribution 
of a fluctuations in a differentially rotating atmosphere. The important difference in 
the present work is that the a fluctuations co nsidered here are statistica l ly stationary 
and h o mogeneous. The numeric al simulations of Brandenburg et ah! ( 20081) : Yousef et al 


(|2008> ); ISingh fc Jingadd (|2015l) demonstrated large-scale dynamo action in a shear flow 
with turbulence that is, on average, non-helical. This problem of the shear dynamo 


















































































































Moffatt drift dynamo 

was taken up in a number of analytical works where the quan tity a was assumed to be 


strictly zero, i.e., 11 vamsnes pomr-wise ooxn m space ana rime uvieeorm nogacnevsKii 
20081: Rogachevskii & Kleeorinl 20081 Sridhar Sz, Subramanian 2009allbl: Sridhar Sz, Singh 

20101; Singh & Sridhar 20111). Based on results of Sridhar & Subramanian (2009a.bl); 

Sridhar & Singh (2010h: ISingh Sz Sridhar 

2011). it was 

realized that the mean 
if it vanishes everywhere 

magnetic field cannot grow if a vanishes strictly (i.e., 

instantaneously). Heinemann. McWilliams Sz Schekochihin ( 

2011): McWilliams ( 

20 ljj; 


Mitra Sz Brandenburg (120121) considered temporal a fluctuations and reported growth of 
second moment of mean magnetic field, i.e., mean magnetic energy. However, it must be 
no ted that the possi bility o f the g rowth of also the first moment existed in the calculations 


of iMitra fe Brandenbura (120121) , where the growth occ urs by the process of negative 


diffusion. This issue was also clarified in a recent work bv ISauire Sz Bhattachariea ( 2015 ) 


who proposed an interesting magnetic shear-current effect as being responsible for shear 
dynamos. 

It was realized in Sridhar fc Sinahl ( 2014 ). and also shown in the present work, that 


the negative diffusion remains the only possibility for driving mean-field dynamos in 
such fluctuating a calculations, so long as these fluctuations are delta-correlated in time. 
On the other hand, the process of the negative diffusion is, in a sense, self limiting, as 
it would increasingly create smaller scale structures, where the necessary assumption of 
the scale-separation cannot continue to be valid indefinitely. This would eventually lead 
to breakdown of the two-scale framework. Given such limitations posed by the process 
of negative diffusion, or in other words, strong a-fluctuation dominated dynamos, it is 
desirable to seek possib il ities of mean -field dynamo action when there are only weak 


a fluctuations. Sridhar & Singh (|2014f) found the possibility of the growth of the mean 


magnetic field also in case of weak a fluctuations when they considered memory effects 
into account. Although many of the previous works included shear in their studies, this 
highlights the importance of fluc t uation s in a , which have indeed been measured in 


simulations of iBrandcnburg et al. ( 2008 ). Estimations of such a fluctuations in various 


astrophysical sources will be of immense value. 

The aim of the present paper is to explore large-scale dynamo action, arising solely 
due to an a that is varying stochastically in space and time, with zero mean. We define 
our model in Section 2. Using the first-order-smoothing-approximation (FOSA), we 
derive an integro-differential equation governing the evolution of the large-scale magnetic 
field in Section 3. This equation is non-perturbative in the a-correlation time r a . We 
first consider the case of white-noise a fluctuations (i.e., T a = 0) without any further 
approximation. Assuming small but non-zero a-correlation time, r a ^ 0, and slowly 
varying mean magnetic field, we simplify the integro-differential equation to a partial 
differential equation (PDE) in Section 4. Without loss of generality we explore in Section 5 
one-dimensional propagating modes and solve the dispersion relation to obtain explicit 
expression for the growth rate function. In Section 6 we study dynamo action due to 
Kraichnan diffusivity and Moffatt drift. We conclude in Section 7. 


2. Definition of the model 

Let us consider a fixed Cartesian coordinate system with unit vectors (ei, e?, ef), where 
x = {x\,X 2 ,x^) denotes the position vector, and t is the time variable. We model small 
scale turbulence in the absence of mean motions as random velocity fields, {v(x,t)}, 
and consider an ensemble with its members corresponding to different realizations of the 
velocity field v. Denoting by ( ) the ensemble av erage which obeys standard Reynolds 
rules (see e.g., Brandenburg Sz Subramanianll2005l) . and assuming that the ensemble has 
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(i) zero mean isotropic velocity fluctuations, (ii) uniform and constant ensemble-averaged 
kinetic energy density per unit mass, and (iii) slow helicity fluctuations, we write: 


(vi) = 0 ; (viVj) = SijV o ; 



f-inj h(x, f) ; 


( 2 . 1 ) 


where v 2 = (v 2 /?>) = two-thirds of the ensemble-averaged kinetic energy density per 
unit mass, and /i(a;,t) = ( v• (Vxd)) /6 = one-sixth of the ensemble-averaged helicity 
density. Let be the size of the largest eddies and r c be the velocity correlation time. 
By slow helicity fluctuations we mean that the spatial and temporal scales of variation 
of fi(x,t) are assumed to be much larger than Iq and r c . 

Let B(x,t) be the meso-scale magnetic field, obtained by averaging over the above 
ensemble. This requires a scale separation such that the typical scales of B(x,t ) are 
much lar ger than In- Then the space-time evolution of B(x,t) is given by the follo wing 
equation (Moffatt 1978 : Krause fo R.adler 198(1 [Brandenburg fc Subramanianll2005ll : 


OB 

— = Vx [a{x,t)B] + t7 T V 2 S; 


V-B = 0. 


( 2 . 2 ) 


where 


a = — 2 T c n(x,t), r]T = g + rj t = total diffusivity , 

77 = microscopic diffusivity, r) t = t c Vq = turbulent diffusivity. (2.3) 


The above simple expressions for the turbulent transport coefficients a and r)t are 
valid only in specialized conditions — specifically, under the first order smoothing 
approximation (FOSA) assuming isotropic turbulence, when the high-conductivity 
limit is considered in conjunction with the so-called “short-sudden” approximation, 
where the velocity correlation tim e r c is much smaller than its turnover time tq 
( Courvoisier. Hughes fc Tobias! 2006h. 

In the “double-averaging scheme”0 being employed here, the helicity fluctuations are 
modelled by fluctuating a, which makes equation (El a stochastic partial differential 
equation. As we are finally interested in scales much larger than the scales of the meso- 
scale field, with quantity a being smooth around meso scale but fluctuating at larger 
scales, we can repeat the averaging procedure to obtain the evolution of large scale field. 
Here, the important step is to consider a superensemble over which a{x, t ) is a statistically 
stationary, homogeneous, random function of x and t, with zero mean, a(x,t) = 0. The 
two-point space-time correlation function of fluctuating a: 


a(x,t)a(x',t') = 2 A(x — x')V(t — t'), with 


poo 

! / V(t)dt = 1 , A( 0 ) = r] a ^ 0 . 

Jo 


(2.4) 


Here rj a is the a-diffusivity, introduced first bv iKraichnanl ( 1976 1. Let us define the 
correlation time for the a fluctuations as, 


= 2 


d 11 V(t). 


(2.5) 


f The concept of double-averagin g has been discussed in a number of previous works (see 
e.g.. IKraichnanl Il976l : ISokolovl 1 19971 1. We refer the reader to Section 11 of iMoffattl (1198311 for 
an excellent account of successive averaging scheme over a number of widely separated scales, 
leading to a successive renormalization of turbulent transport. 



































5 


Mojfatt drift dynamo 

The meso-scale field is split as B = B + b, where B is the large-scale magnetic field 
which is equal to the superensemble-average of the meso-scale field, and b is the part of 
the magnetic field that fluctuates on the meso-scales, simply referred to as the fluctuating 
magnetic field from here onwards. Applying Reynolds averaging to the equation (TOD and 
assuming a(x,t) = 0, we obtain the following equation governing the dynamics of the 
large-scale magnetic field: 

_ _ 

— = Vx£ + rj T V 2 B, VB = 0, (2.6) 


where £ = a(x,t)b(x,t) . (2.7) 

To calculate £, the mean electromotive force (EMF), we need to solve for the fluctuating 
field, b(x,t), whose evolution is determined by: 

db 

— = Vx [aB] + Vx [ab ~ ab\ + 77tV 2 B, 

V-b = 0, with initial condition b(a;,0) = 0. (2-8) 

To keep the analysis simple while providing nontrivial extension to the existing models, 
some simplifying assumptions were made, and therefore it is useful to recall the basic 
limitations of our model. Equation (12.211 assumes a local and instantaneous relation 
between the meso-scale EMF and the corresponding magnetic field. Another limitation 
is the choice of isotropic transport coefficients a and rj t , and it is desirable to understand 
the effects of fluctuations in all components of more general tensorial and . This is 
beyond the scope of the present investigation and will be the subject of a future study. 


3. Equation for the large-scale magnetic field 

To derive a closed equation for the large-scale magnetic field, we first solve for the 
small-scale magnetic field. Following the standard closure technique known as the FOSA, 
or in other words, a quasi linear approach, where we ignore the mode coupling term in 
the fluctuating field equation, we drop the term VX [ab — ab\ from equation (12.81) , but 
retain the r]T'V 2 b term[j] Thus the small-scale magnetic field evolves as: 

b = VxM, V- b = 0, b(x, 0) = 0, (3.1) 



where M(x,t) = a(x,t)B(x,t ) is a stochastic source field. Let us define the spatial 
Fourier transform of, say, b(x,t), denoted by b(k,t ), and its inverse transform as: 


b(fc, t) 



x exp (—i k • x) b(x , t) 


and b(x,t ) 


d 3 fc 

( 2^)3 


exp (i k‘ x) b(k, t ) 

(3.2) 


f SS14 dropped the tjrV 2 !) term too, for simplicity, from the evolution equation for the 
fluctuating magnetic field, while studying the shear dynamo problem. They pointed out that this 
would result in overestimation of growth/decay rates for large wavenumbers. This is confirmed 
later in this work. 
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Fourier transforming equation m, we get after some algebra: 


at 


+ grk 2 ) b = \kxM , k-b = 0, 6(fc,0) = 0. 


where M(k,t) = 


(2tt) 3 


d 3 k'a*(k',t)B{k + k',t). 


(3.3) 


Solution to the equation (13.31) satisfying the constraints, k- b = 0 and b(k , 0) = 0, may 
be obtained by direct integration, which gives 


b(k,t ) = / At' exp [— r]Tk 2 (t — t')] i kxM(k,t') 


We use equation (13.dll to calculate the Fourier transform of the mean EMF: 


(3.4) 


S(k,t) = / d d xexp (—i k- x) £(x,t) = / d 3 a:exp(—i k-x) a(x,t)b(x,t) 


d 3 k' d 3 k" S(k' + k" - k) a{k', t) b{k ", t) 


(2tt) 3 j 

l —r [ d 3 k' d 3 k" 5{k' + k" -k) 
7T Y J 


(2 


x / At' exp [— r/Tk" 2 (t — f')] i k" Xa(k\t)M(k" ,t') . (3.5) 

Jo L - 1 


Above expression for £(k,t) is given in terms of the quantity a(k',t)M(k",t r ), which 
has to be calculated. Using equation (13.31) for M , 


a(k',t)M(k",t') = 


(2tt) 3 


d 3 fc"'5(fc',t)5*(fc ,f') B{k" + k"',t') (3.6) 


is a convolution of the large-scale magnetic field and the Fourier-space 2-point correlator 
of the stochastic a. Using equation m we write the following expression for 2-point 
correlator in Fourier space: 

a(k', t)a*(k'", t') = J d 3 x'd 3 x'" exp (—i k'• x' + i k!"• x'") a(x', t)a(x"', t') 

= 2 T>(t—t') J d 3 x'd 3 x"' exp [— i(fc'- x' — k"'- x'")] A(x' — x'”) . 


Using new integration variables, r = x' — x'" and r' = \(x' + x '"), we get 


a(k', t)a*(k"\ t') = 2V(t — t') J d 3 rd 3 r' exp 

= 2 V(t - t'){ 2tt) 3 8{k' - k"')A{k') 
where A(k) = / d 3 r exp (—ifc- r) A(r). 


-i{k'- k'")-r' -~{k' + k'")-r 


A(r) 

(3.7) 

(3.8) 
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A{k) is the complex spatial power spectrum of a fluctuations, with A(—k) = A*(k) 
because A(r) is a real function. From equations (13.61) and (13.71) we write 

a(k',t)M{k",t') = 2 V(t - t')A(k')B(k'+ k",t'). (3.9) 

When equation (13.91) is substituted in (13.51) we obtain a compact expression for the EMF: 

i(M) = 2jf dsV(s) {u(k,s)xB(k,t- s)} , (3.10) 

~ t' d 3 k' ~ 

where U(k,s)= / -—— exp [—ryrik — k') 2 s] i(fe — k')A(k '). (3.11) 

J (2?r) J 

Fourier transforming equation (12.61) . the equation governing the large-scale field is: 

ffp? 

= \kx£ — r]Tk 2 B , k-B = 0. (3-12) 

Thus the set of equations (13.101) (13.121) describe the evolution of the large-scale magnetic 

held, B. Solving these in full generality is beyond the scope of the present investigation, 
and we study this system of closed equations analytically in useful approximations. 


3.1. White-noise a fluctuations 

Before studying finite memory effects, we consider the exactly solvable limit of white- 
noise (i.e., delta-correlated-in-time) a fluctuations, for which, the normalized corre¬ 
lation function is X\vN(t) = S(t ), the Dirac delta-function. This gives r a = 0 from 
equation (12.51) . implying that the memory effects are ignored. The focus of this section 
is to understand the dynamo behaviour due to white-noise a fluctuations. From equa¬ 
tion (13.101) the mean EMF for white-noise is 

£ wn (M) = U{k,0)xB{k,t ), (3.13) 

which, as expected, depends on the large-scale held at the present time only. The quantity 
U{k, 0) may be simplified as: 


C7(fc,0) = i k 


d 3 k' 

(2tt) 3 


A(k') 



i k'A(k') 


= ifcM(O) — 


&A(0 


3* 


= ikria 




M : 


(3.14) 


£=o 


Here, rj n is the q-diffusivity defined in equation (El, and Vm is the Moffatt drift velocity 
( MofFattl 1 19 78l : Sridhar fe Singh 2014 ). Both, r] a and Vm, are constants by definition. 
Substituting (13.131) in (13.121) and using equation (13.141) . we get the following partial 
differential equation for the large-scale magnetic held: 


dB 

~dt 


[ Ik k 2 + i k-V M ] B = 0 , 


k-B = 0. 


(3.15) 
























Here 
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rjx = TjT — r] a = Kraichnan difFusivity, 




f 

Jo 


a(x, t)Va(x, 0) df = Moffatt drift velocity (3.16) 


are the two constants that determine the behaviour of the large-scale magnetic field. 
Note that the a difFusivity contributes a decrement to the difFusivity , and hence aids 
dynamo action ( Kraichnanl 1976 ; Moffatt 1978 ; Sridhar fe Singhl 12014 ). The solution to 
equation (13.151) is given by 


B(k,t) = Q(k,t) B(k, 0). 


k-B = 0. 


where Q(k, f) = exp j— r]Kh 2 t — i (VM-fc)fj. 


(3.17) 

(3.18) 


Equations (13.171) and (13.181) provide complete solution to the problem of white-noise a 

fluctuations, where the growth or decay of mean magnetic field B is determined by the 
Green’s function Q(k,t). We note some general properties: 

1 . Weak a fluctuations have rj a < tjt so that rjx > 0 . In this case modes of all 
wavenumbers k decay. 

2. Strong a fluctuations have r/ a > t/t so that rjx < 0. This belongs to the 
case when a-diffusivity compensates and overcomes the total (turbulent+microscopic) 
difFusivity rjT ■ This leads to the growth of mod es of a ll wav enum bers k by the process of 
negative diffusion which was first obtained bv iKraichnan ( 1976il . It may be noted that 
the process of the negative diffusion is, in a sense, self limiting, as it would increasingly 
create smaller scale structures, where the necessary assumption of the scale-separation 
cannot continue to be valid indefinitely. This would eventually lead to breakdown of the 
two-scale framework. 

3. The Moffatt drift velocity Vm contributes only to the phase and does not determine 
the growth or decay of the large-scale magnetic field. 

Therefore, the necessary condition for dynamo action for white-noise a fluctuations is 
that they must be strong, i.e., tjk < 0. It is necessary to consider V(t) 6(t) to explore 
memory effects, which will have r a 0. This is studied in the next section. 


4. Large-scale magnetic fields when r a is small 

Now we extend our analysis to include the effects of finite memory of fluctuating 
a. By assuming small r a , we reduce general integro-differential equation, given by 
equations (13.101) (13.121) . to a PDE governing the evolution of large-scale magnetic field 
which evolves over times much larger than r a . Then we first consider the Kraichnan 
problem with nonzero r a , but without Moffatt drift Vm, which is studied in detail in 
Section 6. 


4.1. Derivation of the governing equation 

We note that the normalized time correlation function, 27(f), has a singular limit: i.e. 
lim 27(f) = 27 W n( f) = S(t). Here we wish to consider non-zero but small r a which implies 

T a —> 0 

that the function 27(f) is significant only for times t ^ r Q and becomes negligible for 
larger times. We simplify the mean EMF given by equation (13.101) . together with (13.111) . 
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by solving the time integral for small r a . Since the limit lim £(k,t) = £wn (k,t) = 

T a —>- 0 

U(k, 0)xB(k, t.) is evidently non singular, we make the ansatz that, for small r a , the 
EMF can be expanded in a power series in r a : 

~ ~ ~( 1 ) ~( 2 ) 

£(k,t) = 5wn (k,t) + £ (k,t) + £ (k,t) + ... (4.1) 

~ ~(n) 

where £wn( fc,i) ~ 0(1) and £ ( k,t ) ~ 0(t for n ^ 1. Below we verify this ansatz 

up to n = 1, for slowly varying magnetic fields. 

We wish to determine £(k,t) which is correct to first order in r Q , for t r Q . Since 
T>(s) is strongly peaked for times s ^ r a and becomes negligible for larger s as mentioned 
above, most of the contribution to the integral in (13.1011 comes only from short times 
0 < s < T a . Hence in (13.101) we can (i) set the upper limit of the s integral to +oo; (ii) 
keep the terms inside the { } in the integrand up to only first order in s. Below we first 

work out U(k, s ) and B(k, t — s) correct upto O(s). 

• U(k , s) to O(s) : Taylor expansion of the function U(k, s) gives, 


U(k, s) 


U(k, 0) + s— 


+ 0{s 2 ). 

s =0 


(4.2) 


Let us rewrite U(k,s), correct upto 
U(k,s) = P(k) + 

where P(fe) = U(k, 0) 


and 


Q(k) 


dU 

ds 


5=0 


The integrand in equation (14.51) can 
Q(k) (see the appendix): 


0(s), as, 

sQ(k) (4.3) 

= ikija + Vm (from equation ( 13 . 141) 1 , (4.4) 

r A'iy _ 

= _ir?T J {k k ' )2{k k ' )A{k,) ' (4 ' 5) 

be expanded to obtain the following expression for 


Q(k) = —r) T k 2 (ikr] a + V M ) - 2 r) T {k-V M )k + ikrjrCx + 2ir] T k-C 2 - VtC 3 , (4.6) 


where r] a and Vm are the constants defined in equations (12.41) and (13.161) . respectively, 
which shows that they depend on the spatial correlation function A and its fir st spatial 
derivative. Both these constants appeared already in Sridhar fc Singbl ( 2014h . But C\ 




(scalar), C 2 (dyad) and C 3 (vector) are three new constants given by (see the appendix), 


^ = [ v 2 ^)] c=0 


& = [VV^)]^ = 0 


c 3 = 


= [V 2 {V^)}]^ = 0 , (4.7) 


which depend on second or third order spatial derivatives of A. Thus, in general, the 
mean EMF will be determined by five constants ( t]k , Vm, Cl, C 2l C 3) which can all be 
explicitly found once the form of spatial correlation function A(£) is chosen. To keep the 
analysis simple, we assume that the function >4(£) is such that its second or higher order 
spatial derivatives are negligible at £ = 0, and therefore we ignore the constants C±, &2 
and C 3 in the present work. Thus we write: 


Q(k) = - r]Tk 2 P{k ) - 2r) T {k • V M )k . 


( 4 . 8 ) 
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• B(k,t — s) to O(s) : Taylor expansion of the function B(k , t — s) gives, 

t — s) = ~S{k,t) - g d-BjM) + _ ( 4 . 9 ) 

at 

where it is assumed that |£?| s\dB/dt\ s 2 \d 2 B/dt 2 \ » s 3 \d 3 B/dt 3 \ etc. In 

order to find B(k,t — s ) which is correct upto O(s), we need dB/dt in equation (14.011 
only upto 0(1). Using equation (14.111 together with (14.1411 and (14.1411 in equation (14.1211 
we find, 

= ifcxfwN(fcd) — grk 2 B = — [i k-P + rjTk 2 ] B , (4.10) 

0 ( 1 ) 

which is substituted in equation 62D to obtain: 

B{k , t-s) = B(k , t) [ 1 + s (ifc • P + vrk 2 ) ] + 0(s 2 ) (4.11) 

Then, from the above analysis, we may write, 

U(k,s)xB(k,t — s) = P(k)xB(k,t) + 


dB 

~dt 


(i k-P + r] T k 2 ) P(k)xB(k,t) + Q{k)xB(k,t ) +0(s 2 ). (4.12) 


Noting that £wn (k,t) = P(k)xB(k,t) and using the definition of Q(k) as given in 
equation (14.81) . we get: 

{ } of eqn. (14.101) = £wn( fe, t) + s |(ife • P) £wn — 2 T lT (k-V M )kxB(k,t)}+0{s 2 ). 

( 4 . 13 ) 

Using equation (14.1411 and the properties of T>(t), given in (12.411 and (I2.5L we solve the 
integral over s in equation (14.1011 to obtain the mean EMF: 

£(M) = £wn(M) + T a |(ifc • P) £wn ^ 2 i) T {k-V M )kxB(k,t)} ( 4 . 14 ) 

accurate to 0(r a ), which verifies the ansatz of equation (14.11) up to n = 1, as claimed. We 
note that equation (14.1411 for the mean EMF is valid only when the large-scale magnetic 

field is slowly varying. To lowest order this condition can be stated as: | B | r Q | dB/dt\. 

Using equation (14.1011 for dB/dt , we see that the sufficient condition for equation (14.1411 
to be valid is that the following two dimensionless quantities be small: 

\vk k 2 T a \ < 1 , \kV M r a \ < 1 . ( 4 - 15 ) 

Substituting (14.1411 in (13.1211 we obtain the following PDE governing the evolution of 
large-scale magnetic field: 

^ - ?? 2 fc 4 r„ + (k-V M ) 2 T a ] + i (k- V M ) [ 2(t?« + Vr)k 2 T a - l]|s 

( 4 . 16 ) 


4.2. The Kraichnan problem with non zero r a 
First we consider the Kraichnan problem and extend it to include finite r Q in order 
to understand the combined effect of the a fluctuations when Vm = 0 , but r/ a > 0 and 
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Figure 1. Growth rate 7/cr as a function of \k/k a \ , when Vm = 0 . Weak ( [t/k > 0) and strong 
{r/K < 0) a fluctuations correspond to bold and dashed curves, respectively. Results of SS14 
are also shown; Dash-dotted and triple-dot-dashed curves correspond to weak and strong a 
fluctuations, respectively. 


r a > 0. SS14 defined a length scale whose corresponding wavenumber, 

k a = (ij 0 t„)' 1/2 > 0. (4-17) 


This is used to define high or low wavenumbers (fc); |fc| > k a are called high wavenumbers, 
and |fc| < k a are called low wavenumbers. When Vm = 0, one of the two conditions in 
( 14 . 151 ) are met trivially, and the other one implies that |fc| must be small enough such 
that \r/Kk 2 T a \ -C 1 . Using ( 14 . 171 ) and setting Vm = 0 in equation ( 14 . 161 ) . we get 


dB 

dt 


-ri K k 2 


1 + * 

Vk 



B, 


with k-B = 0. (4-18) 


The solutions are of the exponential form, B{k , t) = Bo(k) exp (yt), where k • Bo(k) = 
0. Substituting this in equation (14.181) . we get the growth rate as: 


7 = 



when \r]Kk 2 T a \ <C 1. (4-19) 


Following SS14 we normalize the growth rates by the characteristic frequency (a) which 
is defined as 


= \m\kl 


\VT - Va 


a 


> 0 . 


(4.20) 
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Below we consider the behaviour of the growth rate as a function of the wavenumber, 
for weak and strong a fluctuations, and show this in Fig. m where we also show the 
corresponding results of SS14 to illustrate some qualitative differences!]] 

Weak a fluctuations: This has 0 < r] a < r]T, so that r]x is positive. The normalized 
growth rate (y/u) may be expressed as 


7 

<7 



\m\ W) 


when 



(4.21) 


The growth rate is negative definite for finite k and it is a monotonically decreasing 
function of the wavenumber; shown by bold curve in Fig. Q. This is even qualitatively 
different from SS14 where the high wavenumbers always grow for weak a fluctuations 
(dash-dotted curve). This also highlights the fact that the inclusion of the resistive term 
in the fluctuating field equation is a nontrivial extension of the SS14 model. 

Strong a fluctuations: This has 0 < r]T < so that t/k is negative. The normalized 
growth rate ( 7 /( 7 ) may be expressed as 


7 

a 



1 - 


\Vk\ 



when 



(4.22) 


In this case the growth rate can be positive for a range of wavenumbers, before becoming 
negative at larger wavenumbers; shown by bold curve in Fig. ©■ This is qualitatively 
similar to the results of SS14 (triple-dot-dashed curve) in this regime and the differences 
at large wavenumbers arise due to reasons mentioned above. 

Thus the necessary condition for dynamo action is that the a fluctuations must be 
strong, i.e., tjk < 0. Recall that this is, in a sense, similar to the white-noise or the 
original Kraichnan model, but we note the following important difference: in white- 
noise case, the growth rate increases monotonically with wavenumber (k) for strong a 
fluctuations, with largest allowed wavenumbers (smallest allowed length scales) growing 
the fastest; whereas, here, we find that the growth rate is a nonmonotonic function of k, 
and as a result, there exists a wavenumber cutoff at some large k beyond which the growth 
rate turns negative. This makes it a special dynamo as the magnetic power at smallest 
length scales would be suppressed due to the existence of wavenumber cutoff (see dashed 
curve in Fig. ©), thus enabling a bonafide large-scale dynamo, unlike the white-noise 
case where much of the magnetic power lies at the smallest allowed length scales. 

It would be indeed interesting if even weak a fluctuations could lead to large-scale 
dynamo action. In the next section we explore the combined effect of Vm 7 ^ 0 and 
T a ^ 0, and ask the following question: when both r a and Vm are nonzero, could large- 
scale magnetic fields grow even when a fluctuations are weak, i. e. when rjx > 0 ? 


5. Growth rates of modes when r a is non zero 

We consider one-dimensional propagating modes for the general case when all the 
parameters {rjonVM,r a ) can be non zero. Below we derive the dispersion relation and 
study the growth rate function. When the wavevector k = (0,0, k) points along the 
“vertical” (± 63 ) directions, B 3 must be uniform and is of no interest for dynamo action. 

Hence we set B 3 = 0, and take B(k, t) = Bflk, t)e 1 + U 2 (/c, f)e 2 . The equation governing 
the time evolution of this large-scale magnetic field is obtained by setting fci j2 = 0 , 

f We note, however, that the expressions for the growth rates in SS14 were simpler and did 
not explicitly depend on the factor g a /\gK\■ 
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= k and B 3=0 in equation (14.161) : 


^ = | [-rj K k 2 - r]lk 4 r a + (kV M 3) 2 T a ] + i kV M 3 [ 2(?7a + Vr)k 2 T a - l]|s (5.1) 

We note that each component of mean magnetic field evolves independently of other 
components. The nature of this dynamo is therefore different from standard a 2 - or 
acc-dynamos where evolutions of various components are coupled with each other thus 
facilitating cross-coupling dynamo. We seek modal solutions of the form, 


B(k,t) 


Bai(k)ei + £?02(fc)e2 


exp (At), 


(5.2) 


and substitute it in equation (15.11) to obtain the following dispersion relation: 

A = {-r) K k 2 - rfck^Ta + (fcI 4 r 3 ) 2 T Q ] + \kV M 3 [2(77 a + r) T )k 2 T a - l] ( 5 . 3 ) 

Of particular interest is the growth rate 7 = Re{A} , because dynamo action corresponds 
to the case when 7 > 0 . From the dispersion relation m we have: 

7 = -r] K k 2 - r] 2 k A T a + ( kV M 3) 2 T a (5.4) 


We refer the reader to Appendix B for some properties of the growth rate function in 
terms of useful dimensionless parameters and to Appendix C for their physical meanings. 


6. Dynamo action due to Kraichnan diffusivity and Moffatt drift 

Now we turn to the most general case when both, Kraichnan diffusivity and Moffatt 
drift, are nonzero, and a-fluctuations have finite correlation times. SS14 defined a new 
time-scale involving tjk and Vm 3 as: 

T* = (\vk\/V 2 3 ) > 0. (6.1) 


We provide below the expressions for dimensional growth rate 7 as a function of the 
wavenumber k , for weak and strong a fluctuations. It turns out that the nature of dynamo 
action depends on whether r a is smaller or larger than t*. 

Weak a fluctuations: This has 0 < rj a < tjt, so that t]k and £k are both positive. 
From equation (15.41) . the dimensional growth rate may be expressed as: 


7 



( 6 . 2 ) 


where the characteristic frequency cr is defined earlier in equation (14.201) . We consider 
following two cases: 

(a) Case r a < t* : In this case the growth rate 7 is negative at all wavenumbers as may 
be seen from the solid curve in the left panel of Fig. ®. 

(b) Case r a > r* : Here the growth rate is positive for a range of wavenumbers and 
it is nonmonotonic function of k. Starting from zero, it first increases with k : attains a 
maximum positive value, 


0 / max 


cr\VK\ \-[a_ 
4?7 a [ T * 


at 


\k\ = k n 




(6.3) 


and then it decreases monotonically for larger fc, turning negative at sufficiently high 
wavenumber; see bold curve in the right panel of Fig. ©. We note that the growth rate 
becomes negative at high enough wavenumbers, thus exhibiting a high wavenumber cutoff, 
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\k/K\ WK\ 

Figure 2. Growth rate 7/cr plotted as a function of \k/k a \ for weak a fluctuations', shown by 
bold curves. Panels (a) and (b) correspond to the case when t<*/t* =0.1 and 10.0, respectively. 
Dashed curves correspond to the results of SS14. 


which would enable a bonafide large-scale dynamo with suppression of magnetic power at 
smaller scales. 

Strong a fluctuations: This has 0 < t]t < rj a , so that r/x and £k are both negative. 
From equation (15.41) . the dimensional growth rate may be expressed as: 


7 = a ■ 


— + 1 

T* 



(6.4) 


As shown by solid curve in Fig. ([3]), the growth rate 7 starts from zero at |fc| = 0, 
increases with k to reach a maximum positive value, 

-12 r 1 1 / N -| 1/2 


/max — 


v\vk\ 

4?7a 


+1 


at 


|^| — ^max — 


\m\ 

‘2f] a 


+ 1 


T* 


(6.5) 


beyond which it begins to decrease monotonically, and becomes negative for sufficiently 
large wavenumbers. 

In_both Figs. © and © we also show the corresponding results of lSridhar fe Singh! 
( 2014h to illustrate some qualitative differences^ We n otice good agreement a t low 
wavenumbers whereas at large wavenumbers the theory of lSridhar fc Singh ( 2014 1 over¬ 
predicts the growth rates. They had already pointed out that such overestimation of 
growth/decay rates at large wavenumbers would be expected as they dropped the 
diffusion term from the evolution equation for the fluctuating magnetic field. In the 
present analysis where we retain this term, we find that it affects the growth rate 7 
in such a way that it now always exhibits a large wavenumber cutoff beyond which 7 
becomes negative. At large wavenumbers, the new predictions for the growth rates are 
even qualitatively different from the results of SS14 model and have only been possible due 
to a nontrivial extension of the previous work. Particularly interesting is the possibility of 
dynamo action in case of weak a fluctuations due to finite Moffatt drift, where a window 
of small to intermediate wavenumbers allows dynamo growth; see right panel of Fig. ©). 


f We note, however, that the expressions for the growth rates in SS14 were simpler and did 
not explicitly depend on the factor g a /\gK\■ 
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Figure 3. Growth rate 7 ja plotted as a function of |fc/fc a | for strong a fluctuations; shown by 
bold curve. Dashed curve correspond to the results of SS14. 


7. Conclusions 

We have developed a theory of large-scale dynamo action where the mean magnetic 
field grows solely due to an a parameter that is varying stochastically in space and 
time with zero mean. Using the first-order-smoothing-approximation or the quasi linear 
approach, we derived a closed integro-differential equation governing the evolution of the 
large-scale magnetic field, which is non-perturbative in the a-correlation time (r a ). This 
is the main result of this paper where we have generalized the Kraichnan-Moffatt (KM) 
model ( Kraichnanl 197 fit Moffatt|ll978 ), to includ e effects of nonzero a-co rrelation time, in 
a spiri t similar to that of Sridhar fc Singhl ( 2014h . We, however, note that Sridhar fc Singhl 
(12014 1 included shear in their analysis while ignoring the diffusion term from fluctuating 
field equation, whereas here we ignore shear but include the diffusion term, which is 
necessary for making comparisons with results from future numerical experiments. We 
show that statistically anisotropic a fluctuations give rise to a drift velocity, called Moffatt 
drift, which contributes a new term in the mean EMF. We first applied our model to 
the exactly solvable case of white-noise (r Q = 0) a fluctuations, in which case the mean 
EMF is identical to the KM model and the evolution of the mean magnetic field depends 
on two constants, namely, t he Kraichnan di f fusivity (vk) and the Moffatt drift ( Vm)- 
We confirm earlier findings (Kraichnan 19761 : lMoffattlfl978 : Sridhar &; Singh 2014li that 
when r a = 0, (a) the necessary condition for dynamo action is that the fluctuations must 
be strong, and (b) the Moffatt drift contributes only to the phase and does not determine 
the growth or decay of the large-scale magnetic field. 

In order to explore memory effects of fluctuating a on dynamo action, we considered 
nonzero r a . Assuming that the r a be small and the large-scale magnetic field is slowly 
varying, we reduce the general integro-differential equation to a partial differential 
equation and state sufficient conditions for its validity. Here each component of the mean 
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magnetic field evolves independently of the other components. We provided an explicit 
expression for the growth rate of the mean magnetic field and studied its behaviour as 
a function of wavenumber (k), for different choices of parameters involved. Some salient 
results may be stated as follows: 

(i) In the absence of the Moffatt drift, the necessary condition for dynamo action is 
that the a fluctuations must be strong. This is, in a sense, qualitatively similar to the 
white-noise or the original Kraichnan model, except that, here, we find the growth rate 
to be a nonmonotonic function of fc, thus exhibiting a wavenumber cutoff beyond which 
it becomes progressively negative. 

(ii) For nonzero r a , the Moffatt drift contributes positively to the dynamo growth and 
it can always facilitate large-scale dynamo action if sufficiently large, even in case of 
weak a fluctuations. 

(iii) In the most general case when both, the Kraichnan diffusivity and the Moffatt 
drift, are nonzero, and r a is finite, we find the possibility of dynamo growth in both 
regimes (weak and strong) of a fluctuations. We also determine the growth rate and 
corresponding wavenumber of the fastest growing mode. 

(iv) We find that there always exists a wavenumber cutoff at some large k beyond 
which the growth rate turns negative, irrespective of weak or strong a fluctuations. This 
makes it a special dynamo as the magnetic power at smallest length scales would be 
suppressed, thus enabling a bonafide large-scale dynamo. 

Thus a minimal extension of KM model to include effects of finite memory results in a 
large-scale dynamo, driven by the Moffatt drift which arises in presen ce of statistical ly 
anisot ropic a fluctuations. Such a possibility was first discussed in Sridhar & Singh 
1 2014( 1. It is particularly intriguing to find that even weak a fluctuations could lead to the 
growth of mean magnetic field due to finite Moffatt drift, where the maximum growth 
occurs at intermediate length scales (~ few k a ). Due to k 2 contribution to the growth 
rate, the Moffatt drift driven dynamo appears to be of negatively diffusive type, with 
coefficient of (negative) turbulent diffusion being VjQ :i T a . However, it is different from 
the usual picture of negative diffusion where the maximum growth occurs at smallest 
length scales. Therefore, while the usual negative diffusion cannot continue indefinitely 
in the mean-field framework for reasons stated earlier, the Moffatt drift driven dynamo 
does not have such limitations, resulting in a bonafide large-scale dynamo action. This 
analysis leading to new contributions to the mean EMF is expected to find applications 
in context to astrophysical dynamos, such as, the disk dynamos, the solar dynamo etc. 
Numerical as well as analytical explorations of this new class of large-scale dynamos, by 
also considering fluctuations in all components of tensorial transport coefficients and 
rjij , will be the focus of a future investigation. 
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Appendix A. Derivation of Equations (14.61) and (14.71) 

To simplify the integral in equation (14.51) . let us first expand (fc — k') 2 (k — k') in the 
integrand and write Q(k) as, 

Q(fe) = -i w |t 2 fc J f^A{k')-P j -2k J 0^ {k .k')A(k') 

I II III 


+k I S^k ,2 A(k') + 2 [ ^-(k.k')k'A(k') - [ ^-k' 2 k'A(k') 1(A1) 


(2t r) 3 


(2tt) 3 


(27r) 3 


iv v 

Writing the spatial correlation function as, 

r d 3 fc' 


VI 


AO = 


(2t r) 3 


exp (i k!• £) Aik 1 ) 


we can express the integrals I-VI in the equation (1 A II) as following: 


I = M( 0 ) ; 


II = -i 


VA(0 


|=o 


III = -i k ■ 


VA (0 


j£=o 


(A 2) 


IV = - 


V 2 M(|) 


H=o ’ 


V = - 


(fe-V)VA(O 


€=o ’ 


VI = i 


V 2 {VA(0} „ (A3) 


J^=d 




Using (E2) in D- we obtain the equation (EH), with definitions of C\ (scalar), C 2 
(dyad) and C 3 (vector) as given in the equation (14.71) . Also, recall that M( 0 ) = r] a and 


VACO 


j£=o 


= Vm- Whereas our model does not specify the form of A(£) which is 


by construct a large-scale quantity, but otherwise an arbitrary function, we restrict our 
present analysis to the limit ffy > £- g-, where l a and fyj- are typical scales of variation 
associated with A and large-scale magnetic field, respectively. In this case the terms 
IV-VI involving C’s in the above expressions can be safely ignored in comparison to the 
rest of the terms in equation CO- 


Appendix B. Dimensionless growth rate function 

The expression for the growth rate as given in equation EH may be written in a 
dimensionless form using the parameters first defined in SS14: 

r = 7 r a , P = r] a k 2 T a , e K = rj K k 2 T a , e M = kV M 3 T a ■ (B 1) 

There is just one constraint involving /3 and Ek , coming from (3 + Ek = VTk 2 T a > 0. 
Thus the parameter ranges are given by, 

0 < /3 < 00, /3 + e k > 0, \e k \ < 1 , |em| < 1 , (B2) 

where the latter two conditions come from equation (14.151) . Multiplying equation (15.41) 
by T a > 0 , we obtain the dimensionless growth rates, 

r = — (ek + / 3 2 ) + (B 3) 

as a function of the 3 dimensionless parameters (/? , Ek , £m). Of these, the two parameters 
{ek , £m) can be taken to be independently specified, taking positive and negative values, 
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£ m P 

Figure 4. Growth rate function F plotted as a function of, (i) Em, for (5 = 0.1 (left panel), 
and (ii) /3, for \em\ = 0.3 (right panel); ek = 0.05 in both the cases. 


so long as their magnitudes are small. But /3 ^ 0 is subject to the constraint (3 + £k > 0 . 
Therefore we can rewrite the conditions of dm) as: 

\ek\ *C 1, \em\ <C 1, 

For £k < 0, have \ek\ < P < oo; For ek > 0, have 0 < /? < oo.(B4) 

The dynamo condition is determined by a surface in three-dimensional parameter space 
(spanned by /3 ,Ek ,e:m ) at which T = 0, separating the dynamo region (with r > 0) 
from the non-dynamo region (with r < 0). In Fig. (U) we plot the growth rate function 
r as a function of single parameter, keeping the other two parameters fixed; left and 
right panels show r as functions of Em and /3, respectively, with positive Ek (= 0.05), 
which corresponds to weak a fluctuations. At fixed Ek and /3, r increases quadratically 
with Em , whereas it decreases quadratically with /? at fixed Ek and Em- We note that the 
Moffatt drift, together with finite correlation time of a-fluctuations (parameterized by 
£m), contributes positively to the dynamo growth and can always facilitate large-scale 
dynamo action if sufficiently strong, even in case of weak a fluctuations. 


Appendix C. Possible physical meanings of the parameters 

If the spatial correlation function Apr) varies over scales of order K_ 4 , then the Moffatt 
drift speed Vm ~ giving from equation (16. Ill , r* ~ |/ — 1| t\/r] a , with factor 

/ = r]T/r] a . Note that r* is a parameter which can be uniquely determined once the 
form of A (r) is specified, whereas r a is defined independently by temporal correlation 
funtion T>(t) using equation (12.51) . Two interesting limits can be sought: (i) when t/t "C r) a , 
t* ~ and (ii) when tjt r] a , t* ~ i-e., for weak a fluctuations, it is larger 

by factor / (with / 1). The wavenumber k a as defined in equation (14.171) signifies 

the inverse diffusion length due to a-diffusivity r] a . Similarly the modified turbulent 
diffusion, or Kraichnan diffusivity rjK has an associated resistive scale with corresponding 
wavenumber, say, kx = 1/,Jr\KT a . The dimensionless parameters j3 and Ek as defined in 
equation m characterise the wavenumbers of the modal mean-field solutions in units 
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of k a and kx, respectively, whereas £m normalizes it by the distance traversed by Moffatt 
drift speed in time r a . 
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